Rapid multiple protein sequence search by parallel and heterogeneous computation.

MOTIVATION: Protein sequence database search and multiple sequence alignment generation is a fundamental task in many bioinformatics analyses. As the data volume of sequences continues to grow rapidly, there is an increasing need for efficient and scalable multiple sequence query algorithms for super-large databases without expensive time and computational costs. RESULTS: We introduce Chorus, a novel protein sequence query system that leverages parallel model and heterogeneous computation architecture to enable users to query thousands of protein sequences concurrently against large protein databases on a desktop workstation. Chorus achieves over 100× speedup over BLASTP without sacrificing sensitivity. We demonstrate the utility of Chorus through a case study of analyzing a ∼1.5-TB large-scale metagenomic datasets for novel CRISPR-Cas protein discovery within 30 min. AVAILABILITY AND IMPLEMENTATION: Chorus is open-source and its code repository is available at https://github.com/Bio-Acc/Chorus.

Improving AlphaFold predicted contacts in alpha-helical transmembrane proteins structures using structural features.

Background: Residue contacts maps offer a 2-d reduced representation of 3-d protein structures and constitute a structural constraint and scaffold in structural modeling. In addition, contact maps are also an effective tool in identifying interhelical binding sites and drawing insights about protein function. While most works predict contact maps using features derived from sequences, we believe information from known structures can be leveraged for a prediction improvement in unknown structures where decent approximate structures such as ones predicted by AlphaFold2 are available. Results: Alphafold2's predicted structures are found to be quite accurate at inter-helical residue contact prediction task, achieving 83% average precision. We adopt an unconventional approach, using features extracted from atomic structures in the neighborhood of a residue pair and use them to predicting residue contact. We trained on features derived from experimentally determined structures and predicted on features derived from AlphaFold2's predicted structures. Our results demonstrate a remarkable improvement over AlphaFold2 achieving over 91.9% average precision for held-out and over 89.5% average precision in cross validation experiments. Conclusion: Training on features generated from experimentally determined structures, we were able to leverage knowledge from known structures to significantly improve the contacts predicted using AlphaFold2 structures. We demonstrated that using coordinates directly (instead of the proposed features) does not lead to an improvement in contact prediction performance.

Improving Inter-Helix Contact Prediction With Local 2D Topological Information.

Inter-helix contact prediction is to identify residue contact across different helices in α-helical integral membrane proteins. Despite the progress made by various computational methods, contact prediction remains as a challenging task, and there is no method to our knowledge that directly tap into the contact map in an alignment free manner. We build 2D contact models from an independent dataset to capture the topological patterns in the neighborhood of a residue pair depending it is a contact or not, and apply the models to the state-of-art method's predictions to extract the features reflecting 2D inter-helix contact patterns. A secondary classifier is trained on such features. Realizing that the achievable improvement is intrinsically hinged on the quality of original predictions, we devise a mechanism to deal with the issue by introducing, 1) partial discretization of original prediction scores to more effectively leverage useful information 2) fuzzy score to assess the quality of the original prediction to help with selecting the residue pairs where improvement is more achievable. The cross-validation results show that the prediction from our method outperforms other methods including the state-of-the-art method (DeepHelicon) by a notable degree even without using the refinement selection scheme. By applying the refinement selection scheme, our method outperforms the state-of-the-art method significantly in these selected sequences.

Targeting of plasmodesmal proteins requires unconventional signals.

Effective cellular signaling relies on precise spatial localization and dynamic interactions among proteins in specific subcellular compartments or niches, such as cell-to-cell contact sites and junctions. In plants, endogenous and pathogenic proteins gained the ability to target plasmodesmata, membrane-lined cytoplasmic connections, through evolution to regulate or exploit cellular signaling across cell wall boundaries. For example, the receptor-like membrane protein PLASMODESMATA-LOCATED PROTEIN 5 (PDLP5), a potent regulator of plasmodesmal permeability, generates feed-forward or feed-back signals important for plant immunity and root development. However, the molecular features that determine the plasmodesmal association of PDLP5 or other proteins remain largely unknown, and no protein motifs have been identified as plasmodesmal targeting signals. Here, we developed an approach combining custom-built machine-learning algorithms and targeted mutagenesis to examine PDLP5 in Arabidopsis thaliana and Nicotiana benthamiana. We report that PDLP5 and its closely related proteins carry unconventional targeting signals consisting of short stretches of amino acids. PDLP5 contains 2 divergent, tandemly arranged signals, either of which is sufficient for localization and biological function in regulating viral movement through plasmodesmata. Notably, plasmodesmal targeting signals exhibit little sequence conservation but are located similarly proximal to the membrane. These features appear to be a common theme in plasmodesmal targeting.

Spheromers reveal robust T cell responses to the Pfizer/BioNTech vaccine and attenuated peripheral CD8+ T cell responses post SARS-CoV-2 infection.

T cells are a critical component of the response to SARS-CoV-2, but their kinetics after infection and vaccination are insufficiently understood. Using "spheromer" peptide-MHC multimer reagents, we analyzed healthy subjects receiving two doses of the Pfizer/BioNTech BNT162b2 vaccine. Vaccination resulted in robust spike-specific T cell responses for the dominant CD4+ (HLA-DRB1∗15:01/S191) and CD8+ (HLA-A∗02/S691) T cell epitopes. Antigen-specific CD4+ and CD8+ T cell responses were asynchronous, with the peak CD4+ T cell responses occurring 1 week post the second vaccination (boost), whereas CD8+ T cells peaked 2 weeks later. These peripheral T cell responses were elevated compared with COVID-19 patients. We also found that previous SARS-CoV-2 infection resulted in decreased CD8+ T cell activation and expansion, suggesting that previous infection can influence the T cell response to vaccination.

In situ cell-type-specific cell-surface proteomic profiling in mice.

Cell-surface proteins (CSPs) mediate intercellular communication throughout the lives of multicellular organisms. However, there are no generalizable methods for quantitative CSP profiling in specific cell types in vertebrate tissues. Here, we present in situ cell-surface proteome extraction by extracellular labeling (iPEEL), a proximity labeling method in mice that enables spatiotemporally precise labeling of cell-surface proteomes in a cell-type-specific environment in native tissues for discovery proteomics. Applying iPEEL to developing and mature cerebellar Purkinje cells revealed differential enrichment in CSPs with post-translational protein processing and synaptic functions in the developing and mature cell-surface proteomes, respectively. A proteome-instructed in vivo loss-of-function screen identified a critical, multifaceted role for Armh4 in Purkinje cell dendrite morphogenesis. Armh4 overexpression also disrupts dendrite morphogenesis; this effect requires its conserved cytoplasmic domain and is augmented by disrupting its endocytosis. Our results highlight the utility of CSP profiling in native mammalian tissues for identifying regulators of cell-surface signaling.

Transcription factor Acj6 controls dendrite targeting via a combinatorial cell-surface code.

Transcription factors specify the fate and connectivity of developing neurons. We investigate how a lineage-specific transcription factor, Acj6, controls the precise dendrite targeting of Drosophila olfactory projection neurons (PNs) by regulating the expression of cell-surface proteins. Quantitative cell-surface proteomic profiling of wild-type and acj6 mutant PNs in intact developing brains, and a proteome-informed genetic screen identified PN surface proteins that execute Acj6-regulated wiring decisions. These include canonical cell adhesion molecules and proteins previously not associated with wiring, such as Piezo, whose mechanosensitive ion channel activity is dispensable for its function in PN dendrite targeting. Comprehensive genetic analyses revealed that Acj6 employs unique sets of cell-surface proteins in different PN types for dendrite targeting. Combined expression of Acj6 wiring executors rescued acj6 mutant phenotypes with higher efficacy and breadth than expression of individual executors. Thus, Acj6 controls wiring specificity of different neuron types by specifying distinct combinatorial expression of cell-surface executors.

KIR+CD8+ T cells suppress pathogenic T cells and are active in autoimmune diseases and COVID-19.

In this work, we find that CD8+ T cells expressing inhibitory killer cell immunoglobulin-like receptors (KIRs) are the human equivalent of Ly49+CD8+ regulatory T cells in mice and are increased in the blood and inflamed tissues of patients with a variety of autoimmune diseases. Moreover, these CD8+ T cells efficiently eliminated pathogenic gliadin-specific CD4+ T cells from the leukocytes of celiac disease patients in vitro. We also find elevated levels of KIR+CD8+ T cells, but not CD4+ regulatory T cells, in COVID-19 patients, correlating with disease severity and vasculitis. Selective ablation of Ly49+CD8+ T cells in virus-infected mice led to autoimmunity after infection. Our results indicate that in both species, these regulatory CD8+ T cells act specifically to suppress pathogenic T cells in autoimmune and infectious diseases.

A new algorithm to train hidden Markov models for biological sequences with partial labels.

BACKGROUND: Hidden Markov models (HMM) are a powerful tool for analyzing biological sequences in a wide variety of applications, from profiling functional protein families to identifying functional domains. The standard method used for HMM training is either by maximum likelihood using counting when sequences are labelled or by expectation maximization, such as the Baum-Welch algorithm, when sequences are unlabelled. However, increasingly there are situations where sequences are just partially labelled. In this paper, we designed a new training method based on the Baum-Welch algorithm to train HMMs for situations in which only partial labeling is available for certain biological problems. RESULTS: Compared with a similar method previously reported that is designed for the purpose of active learning in text mining, our method achieves significant improvements in model training, as demonstrated by higher accuracy when the trained models are tested for decoding with both synthetic data and real data. CONCLUSIONS: A novel training method is developed to improve the training of hidden Markov models by utilizing partial labelled data. The method will impact on detecting de novo motifs and signals in biological sequence data. In particular, the method will be deployed in active learning mode to the ongoing research in detecting plasmodesmata targeting signals and assess the performance with validations from wet-lab experiments.

Single-cell transcriptomes of developing and adult olfactory receptor neurons in Drosophila.

Recognition of environmental cues is essential for the survival of all organisms. Transcriptional changes occur to enable the generation and function of the neural circuits underlying sensory perception. To gain insight into these changes, we generated single-cell transcriptomes of Drosophila olfactory- (ORNs), thermo-, and hygro-sensory neurons at an early developmental and adult stage using single-cell and single-nucleus RNA sequencing. We discovered that ORNs maintain expression of the same olfactory receptors across development. Using receptor expression and computational approaches, we matched transcriptomic clusters corresponding to anatomically and physiologically defined neuron types across multiple developmental stages. We found that cell-type-specific transcriptomes partly reflected axon trajectory choices in development and sensory modality in adults. We uncovered stage-specific genes that could regulate the wiring and sensory responses of distinct ORN types. Collectively, our data reveal transcriptomic features of sensory neuron biology and provide a resource for future studies of their development and physiology.

Temporal evolution of single-cell transcriptomes of Drosophila olfactory projection neurons.

Neurons undergo substantial morphological and functional changes during development to form precise synaptic connections and acquire specific physiological properties. What are the underlying transcriptomic bases? Here, we obtained the single-cell transcriptomes of Drosophila olfactory projection neurons (PNs) at four developmental stages. We decoded the identity of 21 transcriptomic clusters corresponding to 20 PN types and developed methods to match transcriptomic clusters representing the same PN type across development. We discovered that PN transcriptomes reflect unique biological processes unfolding at each stage-neurite growth and pruning during metamorphosis at an early pupal stage; peaked transcriptomic diversity during olfactory circuit assembly at mid-pupal stages; and neuronal signaling in adults. At early developmental stages, PN types with adjacent birth order share similar transcriptomes. Together, our work reveals principles of cellular diversity during brain development and provides a resource for future studies of neural development in PNs and other neuronal types.

Human KIR + CD8 + T cells target pathogenic T cells in Celiac disease and are active in autoimmune diseases and COVID-19.

Previous reports show that Ly49 + CD8 + T cells can suppress autoimmunity in mouse models of autoimmune diseases. Here we find a markedly increased frequency of CD8 + T cells expressing inhibitory Killer cell Immunoglobulin like Receptors (KIR), the human equivalent of the Ly49 family, in the blood and inflamed tissues of various autoimmune diseases. Moreover, KIR + CD8 + T cells can efficiently eliminate pathogenic gliadin-specific CD4 + T cells from Celiac disease (CeD) patients' leukocytes in vitro . Furthermore, we observe elevated levels of KIR + CD8 + T cells, but not CD4 + regulatory T cells, in COVID-19 and influenza-infected patients, and this correlates with disease severity and vasculitis in COVID-19. Expanded KIR + CD8 + T cells from these different diseases display shared phenotypes and similar T cell receptor sequences. These results characterize a regulatory CD8 + T cell subset in humans, broadly active in both autoimmune and infectious diseases, which we hypothesize functions to control self-reactive or otherwise pathogenic T cells. ONE-SENTENCE SUMMARY: Here we identified KIR + CD8 + T cells as a regulatory CD8 + T cell subset in humans that suppresses self-reactive or otherwise pathogenic CD4 + T cells.

Nurturing Undergraduate Researchers in Biomedical Sciences.

Undergraduate researchers are the next-generation scientists. Here, we call for more attention from our community to the proper training of undergraduates in biomedical research laboratories. By dissecting common pitfalls, we suggest how to better mentor undergraduates and prepare them for flourishing careers.

Single-Cell Transcriptomes Reveal Diverse Regulatory Strategies for Olfactory Receptor Expression and Axon Targeting.

The regulatory mechanisms by which neurons coordinate their physiology and connectivity are not well understood. The Drosophila olfactory receptor neurons (ORNs) provide an excellent system to investigate this question. Each ORN type expresses a unique olfactory receptor, or a combination thereof, and sends their axons to a stereotyped glomerulus. Using single-cell RNA sequencing, we identified 33 transcriptomic clusters for ORNs and mapped 20 to their glomerular types, demonstrating that transcriptomic clusters correspond well with anatomically and physiologically defined ORN types. Each ORN type expresses hundreds of transcription factors. Transcriptome-instructed genetic analyses revealed that (1) one broadly expressed transcription factor (Acj6) only regulates olfactory receptor expression in one ORN type and only wiring specificity in another type, (2) one type-restricted transcription factor (Forkhead) only regulates receptor expression, and (3) another type-restricted transcription factor (Unplugged) regulates both events. Thus, ORNs utilize diverse strategies and complex regulatory networks to coordinate their physiology and connectivity.

Cell-Surface Proteomic Profiling in the Fly Brain Uncovers Wiring Regulators.

Molecular interactions at the cellular interface mediate organized assembly of single cells into tissues and, thus, govern the development and physiology of multicellular organisms. Here, we developed a cell-type-specific, spatiotemporally resolved approach to profile cell-surface proteomes in intact tissues. Quantitative profiling of cell-surface proteomes of Drosophila olfactory projection neurons (PNs) in pupae and adults revealed global downregulation of wiring molecules and upregulation of synaptic molecules in the transition from developing to mature PNs. A proteome-instructed in vivo screen identified 20 cell-surface molecules regulating neural circuit assembly, many of which belong to evolutionarily conserved protein families not previously linked to neural development. Genetic analysis further revealed that the lipoprotein receptor LRP1 cell-autonomously controls PN dendrite targeting, contributing to the formation of a precise olfactory map. These findings highlight the power of temporally resolved in situ cell-surface proteomic profiling in discovering regulators of brain wiring.

Transsynaptic Fish-lips signaling prevents misconnections between nonsynaptic partner olfactory neurons.

Our understanding of the mechanisms of neural circuit assembly is far from complete. Identification of wiring molecules with novel mechanisms of action will provide insights into how complex and heterogeneous neural circuits assemble during development. In the Drosophila olfactory system, 50 classes of olfactory receptor neurons (ORNs) make precise synaptic connections with 50 classes of partner projection neurons (PNs). Here, we performed an RNA interference screen for cell surface molecules and identified the leucine-rich repeat-containing transmembrane protein known as Fish-lips (Fili) as a novel wiring molecule in the assembly of the Drosophila olfactory circuit. Fili contributes to the precise axon and dendrite targeting of a small subset of ORN and PN classes, respectively. Cell-type-specific expression and genetic analyses suggest that Fili sends a transsynaptic repulsive signal to neurites of nonpartner classes that prevents their targeting to inappropriate glomeruli in the antennal lobe.

Functional divergence of Plexin B structural motifs in distinct steps of Drosophila olfactory circuit assembly.

Plexins exhibit multitudinous, evolutionarily conserved functions in neural development. How Plexins employ their diverse structural motifs in vivo to perform distinct roles is unclear. We previously reported that Plexin B (PlexB) controls multiple steps during the assembly of the Drosophila olfactory circuit (Li et al., 2018b). Here, we systematically mutagenized structural motifs of PlexB and examined the function of these variants in these multiple steps: axon fasciculation, trajectory choice, and synaptic partner selection. We found that the extracellular Sema domain is essential for all three steps, the catalytic site of the intracellular RapGAP is engaged in none, and the intracellular GTPase-binding motifs are essential for trajectory choice and synaptic partner selection, but are dispensable for fasciculation. Moreover, extracellular PlexB cleavage serves as a regulatory mechanism of PlexB signaling. Thus, the divergent roles of PlexB motifs in distinct steps of neural development contribute to its functional versatility in neural circuit assembly.

Stepwise wiring of the Drosophila olfactory map requires specific Plexin B levels.

The precise assembly of a neural circuit involves many consecutive steps. The conflict between a limited number of wiring molecules and the complexity of the neural network impels each molecule to execute multiple functions at different steps. Here, we examined the cell-type specific distribution of endogenous levels of axon guidance receptor Plexin B (PlexB) in the developing antennal lobe, the first olfactory processing center in Drosophila. We found that different classes of olfactory receptor neurons (ORNs) express PlexB at different levels in two wiring steps - axonal trajectory choice and subsequent target selection. In line with its temporally distinct patterns, the proper levels of PlexB control both steps in succession. Genetic interactions further revealed that the effect of high-level PlexB is antagonized by its canonical partner Sema2b. Thus, PlexB plays a multifaceted role in instructing the assembly of the Drosophila olfactory circuit through temporally-regulated expression patterns and expression level-dependent effects.

pH-Responsive Janus Film Constructed with Hydrogen-Bonding Assembly and Dopamine Chemistry.

Poly(acrylic acid) (PAA) was partially grafted with dopamine (PAA-dopa), and then layer-by-layer assembled with poly(vinylpyrrolidone) (PVPON) to prepare hydrogen-bonded (PVPON/PAA-dopa) n film. Polydopamine (PDA) was deposited on (PVPON/PAA-dopa) n film in the presence of oxidant, and hence the whole (PVPON/PAA-dopa) nPDA film was cross-linked. (PVPON/PAA-dopa) nPDA could be utilized as a platform to produce the free-standing Janus film because of the easy detaching process and various chemical reactivity of PDA layer. Ag nanoparticles were formed on (PVPON/PAA-dopa) nPDA film by electroless metallization. 1 H,1 H,2 H,2 H-Perfluorodecanethiol (PFDT) was used to further modify the film through Michael addition. After detaching from the substrate, (PVPON/PAA-dopa)20PDA/Ag/PFDT exhibits reversible swelling-shrinking behavior as the pH value changes. This free-standing film shows Janus character, one side is hydrophobic, whereas the other side is hydrophilic. In addition, the hydrophobic surface exhibits a surface-enhanced Raman scattering effect, whereas the hydrophilic side does not.